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Abstract 



We have investigated D + tt and D* + ir final states and observed the two es- 
tablished L = 1 charmed mesons, the L»i(2420)° with mass 24211^1 MeV/c 2 
and width 20l^ MeV/c 2 and the D|(2460)° with mass 2465 ±3± 3 MeV/c 2 
and width 28^7^6 MeV/c 2 . Properties of these final states, including their 
decay angular distributions and spin-parity assignments, have been studied. 
We identify these two mesons as the ju g ht = 3/2 doublet predicted by HQET. 
We also obtain constraints on Vs/(Ts + Td) as a function of the cosine of the 
relative phase of the two amplitudes in the L>i(2420)° decay. 



* Permanent address: INP, Novosibirsk, Russia 

f Permanent address: University of Hawaii at Manoa 
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I. INTRODUCTION 



Heavy Quark Effective Theory (HQET) predicts the presence of an approximate flavor- 
spin symmetry for hadrons containing one heavy quark (mg ^> Aqcd) HJH- One of the 
outstanding issues in this theory is whether the charm quark is sufficiently heavy for the 
approximations made in the theory to be valid. One testable prediction of HQET is the 
partial wave structure of the decays of the Dj mesons. The Dj mesons consist of one 
charmed quark (Q) and one light quark (q) with relative orbital angular momentum L. 
When L = 1, there are four states with spin-parity J p = + , 1 + , 1 + and 2 + . In the notation 
introduced by the Particle Data Group ||, these states are labeled Dq, D 1 for both 1 + 
states, and D\, respectively. 

Parity and angular momentum conservation place restrictions on the strong decays of 
these D j states to Dn and D*n: the + state can decay only to Dn through S-wave decay, 
either 1 + state can decay to D*n through S-wave or D-wave decays, and the 2 + state can 
decay to both Dn and D*n only through D-wave decays. 

In the decay Dj — > D*n, the helicity angular distribution of the D* can be used to 
analyze the spin of the parent Dj. The helicity angle, denoted by a, is defined as the angle 
between the n~ from the decay D° ; — > D* + ir~ and the n + from the decay D* + — > D°n + , 
both measured in the D* + rest frame. Regardless of the initial polarization of the Dj states, 
the predicted helicity angular distributions are: 



sin 2 a (2 + state) 

? .. v . 1 (pure S-wave 1 + state) (1) 

' (t,sn I i -)_ 3 cos 2 a (pure D-wave 1 + state) 



d N 

oc 



Two Dj states have been observed previously (§-§ . However, the decay angular analyses 
performed were all incomplete and not in total agreement with each other. 

In the limit tuq — > oo, the mesons are described by j = S q + L (total angular momentum 
of the light quark). When mq is large, but not infinite, the properties of the mesons will 
also depend on Sq (spin of the heavy quark) and J = j + Sq (total angular momentum). 
The mesons with L = 1 are then labeled as: 



T U)P _ j Pt> m+ , A (1/2)+ (j = 1/2 doublet) 
J ~ 1 P$"*> + , P 2 (3/2)+ (j = 3/2 doublet) 



The j = 1/2 mesons are predicted to decay exclusively in an S-wave, while the j = 3/2 
mesons decay only in a D-wave |T]|| . Mesons which decay via a D-wave are predicted to be 
relatively narrow (widths of tens of MeV/c 2 ), while those of the other doublet are predicted 
to be quite broad (hundreds of MeV/c 2 ). The p^ 2 ^ state can decay to both Dir and D*ti. 
Models [|9|-|rTH predict the ratio of the branching fractions, BiD^ — ► Dn) / B{D\ —>■ D*ir), to 
lie in the range from 1.5 to 3.0. 

(3/2)4- 

Because the mass of the charmed quark is not infinite, the P| — > D*tt decay can 
also proceed via an S-wave. HQET predicts that this S-wave amplitude is small compared 
with typical S-wave amplitudes, but makes no prediction for the relative magnitudes of the 
S-wave and D-wave amplitudes for this decay. However, particular models do make such 
predictions jlO|JTT . 
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The large data sample collected with the CLEO II detector allows us to identify clearly 
two of the Dj states, to measure their corresponding angular distributions and to test some 
of the HQET predictions. 

II. DATA SAMPLE AND EVENT SELECTION 

The data used in this analysis were selected from hadronic events produced in e + e~ 
annihilations at CESR and collected with the CLEO II detector. Both neutral and charged 
particles are measured with excellent resolution and efficiency by the CLEO II detector. A 
detailed description of the detector can be found elsewhere |L2| . The center-of-mass energies 
used in this analysis were at the mass of the T(4S), Ec.m. = 10.580 GeV, and in the nearby 
continuum. The data corresponds to an integrated luminosity of 1.7 fb _1 . 

We selected events that have a minimum of three charged tracks, a total visible energy 
greater than 15% of the center-of-mass energy (this reduced contamination from two-photon 
interactions and beam-gas events), and a primary vertex within ±5 cm in the z-direction 
and ±2 cm in the r-0 plane of the nominal collision point. All charged tracks were required 
to have dE/dx information. If available, time-of- flight information was also used. 

When reconstructing tt° candidates, we used pairs of photons from the barrel region with 
| cos 6*| < 0.707, where the energy resolution is best. The photons were required to have a 
minimum energy of 50 MeV and to be isolated from charged tracks. All tt° candidates, 
whose invariant mass was within 50 MeV/c 2 of the 7r° mass, were kinematically fit to the 
known tt° mass and were required to have a minimum momentum of 65 MeV/c. 

III. D£° -» D + 7T 

To reconstruct D^ — > D + n~ |, we first reconstructed the D + in the decay mode 
K~7T + 7r + . We required that the decay angle 9k, the angle between the direction of the 
D + momentum in the lab and the direction of the K~ momentum in the D + rest frame, 
satisfy the condition cos 9k < 0.8, since the background peaks near 1. Each D + candidate 
was then combined with each remaining 7r~ in the event. In order to reduce combinato- 
rial background we applied the cuts of x p (Df) = p(Df) / ^/ E^ eam - M(Df) 2 > 0.65 and 
cos^t,- > —0.8, where the decay angle, 9 n , is defined as the angle between the direction of the 
D^ momentum in the lab and the direction of the n~ momentum in the D^ rest frame. We 
then calculated the total probability, P tot , of the candidate using the particle identification 
(dE/dx and time-of-flight) and the reconstructed D + mass. PtotiXtoti -^W) is defined as the 
probability to observe x 2 > Xtot f° r degrees of freedom. The data sample is highly 
contaminated due to the small D + signal to background ratio, and this produces a large 
and broad peaking in the P to t distribution at P tot = 0. For the signal this distribution is 
expected to be flat. We accordingly required P to t > 0.4. 



^References in this paper to a specific state or decay will always imply that the charge-conjugate 
state or decay has been included as well. 
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The spectrum of the mass-difference, M(D + tt~) — M(D + ), for all D + n~ combinations 
surviving the above cuts is shown in Fig. |]. This spectrum was fitted with a third-order 
Chebychev polynomial for the background and a Breit-Wigner resonance shape, convoluted 
with a Gaussian resolution function, for the signal. The a of this Gaussian function was 
fixed to 4.5 MeV/c 2 , as determined from Monte Carlo studies. The region from 380 to 430 
MeV/c 2 was excluded from the fit because this region is populated by feed-down, caused 
by not reconstructing neutrals in the decay chain, Dj — ► D* + n~, with D* + — ► D + tt° or 
£>+7. Our fit yielded 486±ng signal events with a value M(Df) - M(D + ) = 596 ± 3 ± 3 
MeV/c 2 , which corresponds to a D^ mass of 2465 ± 3 ± 3 MeV/c 2 , and an intrinsic width 
T = 28^7^5 MeV/c 2 . The second error is systematic and was estimated by varying the cuts, 
the background parameterization, and the spin of the Breit-Wigner distribution used. Our 
results for the mass and width of this state, along with previous measurements, are listed 
in Table |. 

IV. D<} -> D*+7T 

To reconstruct Dj — > D* + tt~ we first reconstructed D 0, s in the decay modes K~tt + , 
K~7T + 7r°, and K~tt + 7t + tt~ . The D* + candidates were reconstructed by combining each D° 
candidate with each remaining 7r + in the event. Each D* + candidate was then combined 
with each remaining 7r~ in the event. In order to reduce combinatorial background, a cut of 
Xp(Dj) > 0.6 and a cut of cos^ > —0.7 were applied. We calculated P to t using the particle 
identification (dE/dx and time-of-flight), the D° mass, the ir° mass, and the mass-difference 
M(D* + ) - M(D°). A purifying cut of P tot > 0.05 was then imposed. 

The expected angular distributions, outlined previously, can be used to separate the two 
Dj states decaying to D* + ir~. In order to resolve the D\ state and improve the signal to 
background ratio, the Df state was suppressed by requiring | cosa| > 0.8. The M(D* + n~) — 
M(D* + ) spectrum for all combinations passing the above cuts is shown in Fig. ||](a). The 
prominent peak observed is due to the D® state and the shoulder at higher mass is due to 
the Z?2° state- This spectrum was fitted with a fourth-order Chebychev polynomial for the 
background and two Breit-Wigner resonance shapes convoluted with Gaussian resolution 
functions for the signals. The a's of these Gaussian functions were fixed to 4.0 MeV/c 2 , as 
determined from Monte Carlo studies. The mass and width of the higher mass convoluted 
Breit-Wigner were constrained to the measured values obtained above from the decay D^ — > 
D + 7r~ , while the parameters of the other convoluted Breit-Wigner were left free. For the 
suppressed Df state, we obtained 48 ± 30 signal events. For the D\ state, the fit yields 
2861^ signal events, M(D°) - M(D*+) = 411^1| MeV/c 2 , and T = 20±|±1 MeV/c 2 . The 
measured mass-difference corresponds to a D® mass of242ll^MeV/c 2 . The second error is 
systematic and was estimated by varying the cuts, the mass and width of the fixed state, the 
background parameterization, and the spin of the Breit-Wigner shapes used to describe the 
signals. Our results for the mass and width of this state, along with previous measurements, 
are listed in Table [TT|. 

The spectrum of the mass-difference, M(D* + tt~) — M(D* + ), with no cut on the helicity 
angle is shown in Fig. 0(b). A fit to the mass- difference distribution with no helicity angle cut 
yielded M(D^) - M(D* + ) = 412 ± 1 MeV/c 2 and T = 24±l MeV/c 2 in excellent agreement 
with those obtained above. 
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V. FRAGMENTATION 



The momentum spectra of Z^(2460)° -> D+n- and of 0^2420)° - 
obtained by fitting the observed mass-difference distribution in five x p 



Fitting the Peterson fragmentation function [13 



dN 

doc f) 



i 



oc 



x p [l 



-> D* + 7i were also 
bins from 0.5 to 1. 



(3) 



l-Xr, 



to the acceptance-corrected momentum spectra, shown in Fig. Q(a) and (b), we find e p = 
O.034±8:8ii±8:ooi and 0.015±g:88f±g;881, respectively, for the L>*(2460)° and the J D 1 (2420)°. The 
systematic error was estimated by varying the number and size of the x p bins used, and by 
varying the mass and the width of the Dj states. Both spectra are quite hard compared 
with those observed for continuum production of D and D* mesons at these center of mass 
energies [0,0 . 



VI. CROSS SECTIONS AND PRODUCTION RATIOS 

As will be discussed below, measurements of the production rates for the Dj states are, 
in general, sensitive to the Dj alignment, which is uncertain. When the detector acceptance 
is flat in cos 9 n , the dependence on the alignment can be removed either by including events 
from the full range of cos^ or by selecting events with cos^ > 0. Monte Carlo studies 
showed that the efficiency for reconstructing Z?2(2460)° — > D + n~ does not vary over the 
range —1 < cos^ < +1, while that of Dj — > D* + n~ does not vary significantly over the 
range cos 9 n > 0. We accordingly removed the cos9 n cut for the decay of Z^2(2460)° — > D + n~ 
and required cos 9 n > for the decays of Dj — > D* + n~. The measured yields for the decay 
modes D*(2460)° -> D+tt", D*(2460)° -> D*+n-, and 0^2420)° -> D* + n- were then 
determined to be 513 ± 81, 164 ± 41, and 536 ± 43 signal events, respectively. Using the 
fragmentation functions to extrapolate to zero momentum, and after correcting for the 
efficiencies and the relevant D°, D + , and D* + branching ratios, we have extracted the 
production cross sections times the branching ratios for x p > 0: 

v(e+e~ ^ I%X) • B{0>). (4) 

Our measurements, as well as those of ARGUS MM, are summarized in Table [IT]. In 
these calculations we have used the new CLEO II Jl6| measurements of the D* + branching 



ratios. In making comparisons with previous measurements, we scaled the old numbers to 
compensate for the increase of the B(D* + — ► D°7T + ). The systematic errors include the 
uncertainties in the yields of the Dj states, the uncertainties in the branching ratios, and 
the uncertainty in the extrapolation to x p = 0. 
We also extracted the production ratios: 

N(D°j - D+7T') Xp>0 , 6 I N{D+) Xp>0 , 6 (5) 

and 

N(D°j - D*+n-) Xp>0 . 6 I N(D*+) Xp>0 . 6 . (6) 
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Our measurements, as well as those of ARGUS H and CLEO 1.5 H7I, are summarized in 



Table |IV. 

We then determined the ratio of the branching fractions of the D2(2460)° state: 
S[ J D*(2460)° -> D + iT- 



S[D|(2460)° -> D*+n- 



2.2 ±0.7 ±0.6. (7) 



ARGUS H and CLEO 1.5 § have measured this quantity to be 3.7± 1.4± 1.9 and 2.8± 1.0, 
respectively. Our new result agrees well with the HQET predictions discussed previously. 



VII. HELICITY ANGULAR DISTRIBUTIONS 

The various spin-parity hypotheses and the corresponding helicity angle distributions, 
for the general case of decays to D*tc, are listed in Table |V|. 

A study of the helicity angular distribution, cos a, lends support to the identification of 
these states as members of the j = 3/2 doublet. To study the helicity angle distributions of 
J D*(2460)° -> D* + n~ and Z^^O) -> D*+n~, we fitted the M(D* + 7t~) - M(D* + ) mass- 
difference spectra in five bins of cos a from —1 to ±1. The masses and widths of both states 
were fixed to our measured values given above. 

In the decay of the Z^2(2460)°, the D* + and the 7r~ are emitted in a relative D-wave. This 
requires the D* + to have a helicity of ±1. Thus the form of the helicity angular distribution 
for the _D2(2460)° state will be sin 2 a, independent of the initial alignment of this state. 
Monte Carlo studies showed that our efficiency does not vary significantly over the helicity 
angle. The number of D2(2460)° events obtained in each cos a bin is shown in Fig. |](a). 

The general form for the joint decay angular distribution for the Z} 1 (2420) — > D* + n~ is 
given below, Eq. (8). When the detector acceptance is flat in cos 6^, integration over cos#„- 
in either the range — 1 < cos 6^ < ±1, or in the range cos 6^ > 0, will remove dependence of 
the Z} 1 (2420) helicity angular distribution on the alignment of the Z} 1 (2420) state. Monte 
Carlo studies showed that, in the cos 6^ > range, our efficiency does not vary significantly 
over the plane of the helicity angle and the decay angle. The only case where the efficiency 
is a little lower, by ~ 10%, is when cos a is backward and cos^ is forward. In this case, 
the momentum of the slow 7r + in the decay D* + — > D°tt + has its minimum value. We thus 
required cos^ > and corrected for the relative efficiency of the point at cos a < —0.6, and 
included the uncertainty in the alignment of the Z)i(2420)° state in the systematic error of 
the efficiency correction. This yielded the number of Di(2420)° events in each cos a; bin as 
shown in Fig. f|(b). In addition to the statistical errors, the error bars include systematic 
errors due to the uncertainties in the yields and due to the efficiency correction for the point 
at cos a < —0.6. 

We evaluated the \ 2 P er degree of freedom, x 2 /Ndof, and confidence level for various 
hypotheses for the shape of these distributions. The results are listed in Table [VI]. For the 
Z^2(2460)° state the sin 2 a hypothesis is preferred, although the isotropic hypothesis is also 
acceptable. For the D 1 (2420)° state, the cos 2 a hypothesis is excluded for the first time, 
at more than 99% CL. This excludes an alternative interpretation of this state as a radial 
excitation of the D° with J p = 0~ . Because the sin 2 a hypothesis is also excluded, the cos a 
distribution, alone, restricts the Di(2420)° state to the quantum numbers J p = 1 + , 2~, 3 + . . . 
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Since it is difficult to produce L > 1 in the fragmentation process, J p = 1 + is the preferred 
assignment. Because the D 1 (2420)° width is relatively small, this state is identified as the 



(3/2)- 



Pf / J+ and the Z^(2460)° as its doublet partner, the P 

To compare with previous results, we also fit the distribution for the D^^AQQ) decay to 
the form, A(l + B cos 2 a), which gave the result, B = -0.74t£i with x 2 /N do f = 0.6/3 and 
a CL = 90.6%. Fitting the distribution for the Di(2420)° to the same functional form gave 
the result, B = 2.74t£$j with x 2 /N do f = 2.2/3 and a CL = 53.2%. 

Assuming that the Di(2420)° does indeed have J p = 1 + , both S and D wave amplitudes 
are allowed in its decay to D* + tt~. The general form of the _D 1 (2420)° — > D* + tt~ angular 
distribution is then: 



dN 



sin 2 a 



d cos a d cos 9 n 



+ 



cos 2 a 



(1 + cos 2 9 W ) + poo(l - 3 cos 2 ff )] [2T S + T D + 2^/2T s T D cos 

"(l-cos 2 ^)-p o(l-3cos 2 ^)] {V s + 2T D - 2 y/2T s T D cos<^}, (8) 

where a is the helicity angle, 9 n is the decay angle, Ts is the S-wave partial width, is the 
D-wave partial width, if is the relative phase of the two amplitudes, and poo is the fraction 
of Z)i(2420)° with helicity in the lab frame. Integrating over cos 6^ from -1 to +1, or from 
to +1, gives: 



1 dN 
N d cos a 




R) 



1 + 3 cos 2 a 



+ a/2R(1 - R) cosip 



1 — 3 cos a 



(9) 



where R=r 5 /(r5 + T D ). Had we not released the cos^ cut or applied a cos^ > cut, 
the helicity angle distribution would also have depended on the alignment of the initial 
state. Previous analyses |)],[7||| did not remove this cut. Once R and cosy? are specified, the 
shape of the expected cos a distribution is fixed, and one can then determine the \ 2 th a t 
the distribution fits the data points in Fig. f|(b). The shaded region in Fig. || shows the 
region of the R-cosy? plane which is allowed at the 90% CL. The allowed regions fall into 
two categories. First, if R is small then all values of cos<p are allowed. Second, if cos<^ is 
negative, then a large S-wave partial width is allowed. This can result from mixing between 
the two 1 + states and it may explain the difference between the measured value of the ratio 



V(Df)/V{Dl) ~ 1, and the HQET |MT| prediction ~ 3. 



VIII. CONCLUSIONS 

In conclusion, we have observed the two charmed mesons of masses 2465 ± 3 ± 3 and 
2421±|±1 MeV/c 2 , and of widths 28l^ and 20±f±| MeV/c 2 , respectively. The observed 
helicity angular distributions are in good agreement with expectations for the L = 1 cu 
£>£(2460)° and Di(2420)° states. We have made the first measurement of T S /{T S + T D ) as 
a function of the relative phase of S and D wave amplitudes in the Di(2420)° — > D* + tt~ 
decay. In addition, the measured widths of both states are relatively narrow, consistent 
with the predictions for D-wave decays of these states. We have also determined the ratio 
of the branching fractions of the two decay modes of D2(2460)° into D + tt~ and D* + tc~, 
obtaining a value of 2.2 ± 0.7 ± 0.6 which is consistent with the expectation from HQET. 
Taken together, these results constitute strong evidence for identifying these two Dj states 
as the P 2 (3/2)+ and P^ ,2)+ (j = 3/2) doublet predicted by HQET. 
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TABLES 



TABLE I. D|(2460)° mass and width. 



Experiment Mass (MeV/c 2 ) Width (MeV/c 2 

CLEO II 2465 ± 3 ± 3 28+f+^ 

CLEO 1.5 2461 ± 3 ± 1 2 ^-li-io 

ARGUS 2455 ± 3 ± 5 ^-ll-lo 

E691 2459 ±3±2 20 ±10 ±5 

E687 2453 ±3±2 25 ±10 ±5 



TABLE II. L>i(2420)° mass and width. 



Experiment Mass (MeV/c 2 ) Width (MeV/c 2 

CLEO II 2421+^ 20^1+1 

CLEO 1.5 2428 ± 3 ± 2 23±|±J° 

ARGUS 2414 ± 2 ± 5 13±|1J° 

E691 2428 ± 8 ± 5 58 ± 14 ± 10 

E687 2422 ±2±2 15±8±4 



TABLE III. Summary of Dj cross sections times branching ratios: <r(e + e — > DjX) ■ B(Dj). 

Decay Mode CLEO II (pb) ARGUS (pb) 

1^(2460)° -> D+-K- 21.4 ±3.4 ±4.2 68 ± 21 ± 28 

L>|(2460)° -> D*+ir- 9.5 ±2.4 ±1.8 19 ± 4 ± 6 

£>i(2420)° -» D* + vr- 28.5 ±2.3 ±3.6 32 ± 4 ± 9 



TABLE IV. Summary of D°j production ratios: N(D°j -> r> W+ 7r-)z P >0.6 / A r (-D ( * )+ ) a;p 



Decay Mode 


CLEO II (%) 


CLEO 1.5 (%) 


ARGUS (%) 


£>5(2460)° -> D+vr- 


4.7 ±0.7 ±0.7 


in +2+2 

iu -2-l 


1 1 +4+5 


L>^(2460)° D*+ir- 


2.1 ±0.5 ±0.4 


O R + l.0 + 0.4 

°-°-1.0-0.8 




L>i(2420)° -» L>*+vr- 


6.8 ±0.6 ±0.9 


O Q + l.1 + 0.5 

°- a -l. 1-0.5 
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TABLE V. List of spin-parity hypotheses and the corresponding helicity angle distributions. 
A\ is the amplitude to produce D* with helicity A. 

Spin-Parity Hypothesis Angular Distribution 

+ forbidden 

0~ cos 2 a 

l-,2 + ,3",... sin 2 a 

l+,2-,3+,... |^io| 2 sin 2 Q+ U o| 2 cos 2 a 



State 



TABLE VI. x 2 an d CL for various angular distributions. 



Angular Distribution 



X 2 /N dof 



CL 



L>|(2460)° 



£>i(2420)° 



sin 2 a 
isotropic 



1 + 3 cos 2 a 
isotropic 
cos 2 a 
sin 2 a 



1.2/4 
2.5/4 

2.3/4 
23.9/4 
28.2/4 
93.2/4 



87.8% 
64.5% 

68.1% 
83.6 x 1(T 6 

11.4 x 1(T 6 

27.5 x 1(T 20 
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FIGURES 



FIG. 1. Fhe M(D+tt-) - M(D + ) mass-difference distribution. 



FIG. 2. The M(D* + ir-) - M(D* + ) mass-difference distribution for (a) |cosa| > 0.8 and (b) 
— 1 < cos a < +1 . 
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FIG. 3. The momentum spectra of (a) L>?;(2460) -» D+w' and (b) Z) 1 (2420)° -> L>* + vr-. 
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FIG. 4. The normalized helicity angular distributions for (a) Z?2(2460)° decay and (b) 
£>i(2420)° decay. 



FIG. 5. Plot of R = Ts/ (Ts+Td) versus cosine of the relative phase of S and D wave amplitudes 
in the Z) 1 (2420) decay. The shaded area represents the 90% confidence level allowed region. 



15 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403359vl 
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